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ABSTRACT 

We present predictions of centimeter and millimeter radio emission from reverse shocks in the early 
afterglows of gamma-ray bursts with the goal of determining their detectability with current and future 
radio facilities. Using a range of GRB properties, such as peak optical brightness and time, isotropic 
equivalent gamma-ray energy and redshift, we simulate radio light curves in a framework generalized 
for any circumburst medium structure and including a parametrization of the shell thickness regime 
that is more realis tic than the simple a ssum ption of thick- or thin- shell approximations. Building 
on earlier work by Mundell et al. (20071 and Melandri et al. (2010) in which the typical frequency 
of the reverse shock was suggested to lie at radio, rather than optical wavelengths at early times, 
we show that the brightest and most distinct reverse-shock radio signatures are detectable up to 
0.1—1 day after the burst, emphasizing the need for rapid radio follow-up. Detection is easier for 
bursts with later optical peaks, high isotropic energies, lower circumburst medium densities, and at 
observing frequencies that are less prone to synchrotron self-absorption effects - typically above a few 
GHz. Given recent detections of polarized prompt gamma-ray and optical reverse-shock emission, we 
suggest that detection of polarized radio/mm emission will unambiguously confirm the presence of 
low-frequency reverse shocks at early time. 

Subject headings: gamma-ray burst: general — radiation mechanisms: non-thermal — radio continuum: 
general 


1. INTRODUCTION 

In the standard fireball model of gamma-ray bursts 
(GRBs), internal shocks in the expanding flow produce 
the prompt q-ray emission that characterizes a GRB. As 
the expanding ejecta collide with the surrounding cir¬ 
cumburst medium, a fading afterglow is produced, which 
is comprised of two components: a forward shock (FS) 
that propagates outwards into the ambient medium and 
a reverse shock (RS) that travels backwards into the 


on-co ming flow (Meszaros et al. 1994 Meszaros & Rees 


1994). It was expected that bright optical hashes, pro- 
duced by the RS similar to that found in GRB 990123 


(lAkerlof et al. 

1999 

Sari & Piran 

|1999a 

Kobayashi & 

Sari 2000 Ko 

bayas 

til 

2000), woul 

d be common in the 


Despite ten years of accurate GRB localizati ons dis¬ 
seminated a utomatically in real-time from Swift (Gehrels 
et al. 2004) and rapid, ground-based fol low-up by au¬ 
tonomous robotic optical telescopes (e.g.,iMonfard ini et. 
al. 1120061 iGuidorzi et al.||2011| |Virgili et al.||2013(), only 
action (^hVo) of early-time optical light curves 


a small 


show clear evidence of optical RS emission ([Meszaros & 

Rees 

1999[[Sari & Piran|1999bl[Kobayashi|2000t|Rommg 

et al. 

20061 [Gomboc et al. 120081 [Kopac et al. 120131 IJapeli 

et al. 

20l4[). 


Mundell et al. 

(2007 

) suggested that this lack of bright 


chrotron frequency of the RS already lies at radio fre¬ 


quencies at early time. Melandri et al. (2010) extended 
this low frequency scenario to a sample of 19 GRBs with 
well-sampled optical light curves, each with a single op¬ 
tical peak that was consistent with the typical frequency 
of the FS lying close to the optical band. They pro¬ 
duced model radio light curves for FS and RS emission, 
including a simple parametrization of synchrotron self¬ 
absorption at early time. 

Radio interferometers, such as the Jansky Very 
Large Array (JVLA), the Atacama Large Millime¬ 
ter/submillimeter Array (ALMA) and pathfinders for the 
Square Kilometer Array, will come on line or have been 
significantly upgraded to provide unprecedented sensi¬ 
tivity, providing mjy and fiJy level observations in rel¬ 
atively short integ ration times and at ever-shorter re¬ 


sponse time s (e.g., Chandra & Frail 2012 Laskar et al. 

2014). Ultiniately, sensitive radio surveys of the 
transient sky will become routine with future facilities 
such as the Square Kilometer Array (SKA). Estimates 
for radio GRB event rates are beginning to emerge but 
these currently neglect RS emission, focusing instead on 
FS emission, for which radio lig ht curves have tradition¬ 
ally been better characterized (fChirlan da et al. 2013 
~ ,ed by 


et al. (2010 
GRB radio 


Burlon et al. 2015 [Metzger et aT 
this, and fol lowing [Mundell et al. 


^3). Motive 

]2007) and Melandri| 


we provide new theoretical predictions of 


ight curves including RS emission within an 
updated low-frequency model framework, incorporating 
corrections from hydrodynamical simulations and a gen- 
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eralized circumburst medium. The latter addition stems 
from the practice of analyzing GRBs in a constant ISM 
(fc = 0) and wind-type medium (fc = 2), but also from 
modeling of GRBs with structured optical e mission that 
imply an intermedia te stratihcation {k ^ 1, Liang et al. 


2013 


Yi et al.||2013) and detailed hydrod ynamical light 


curve modeling (e.g. de Golle et al. [20121. 


In Section 2 we present the updated low-frequency 
model formalism, in Section 3 we detail our Monte Carlo 
simulations, in Section 4 we present our results, in Sec¬ 
tion 5 we discuss the implications of the simulations, and 
in Section 6 we highlight the conclusions. 

2. LOW-FREQUENCY MODEL FRAMEWORK 

For GRBs with well-sampled early optical afterglows 
we use a set of simple assumptions about the relationship 
between the RS and FS emission at the deceleration time 
to construct predicted lo w-frequency light curves similar 
to Melandri et al. (2010) but generalized for any circum- 


burst medium structure. The subset of bursts with a 
single peak in their optical light curves are of particular 
interest because they provide information about the time 
of deceleration of the fireball, if we assume tp,opt ~ idcJM 
Using the deceleration time estimated from the opticm 
light curve and the high-energy properties of the burst 
we estimate the Lorentz factor Fq for a general stratified 
medium, where the circumburst medium deri sitv is given 


1999 


2000) and where 


by n = AR~^ (jChevalier & Li 
i? is a distan ce irom the prog enitor star. We adopt the 
formalism of Yi et al. ( 2013 ) where n = no{R/Ro)~^, 
with typical values of ng ~ 1 cm“^ and Rq ^ 10 ^^ cm. 
The Lorentz factor is approximated by: 




(3 - k)E^^iso 


47r TyngR^mpC^ 

k ^ ^p,opt / 


-k)l l/(8-2fe) 


( 1 ) 


where c is the speed of light, nip is the proton mass, rj is 
the radiative efficiency, E^ iso is the isotropic equivalent 
energy of prompt gamma-ray emission, and Ct oc 
is a numerical factor that corrects for the shell thick¬ 
ness regime, as d etermined from numerical h ydrody- 
namic simulations (|Harrison & Kobayashi||2013t 2015 in 
prep). 

At the deceleration time, the RS peak frequency and 
peak spectral flux density are estimated by: 


C^rnTg ZZm,f 


( 2 ) 


F, 


i/'.max.r 


Gi^^FoFb jnax.f 5 


( 3 ) 


where k) and k) are numerical factors that 

correct for the shell thickne ss regime for k = [0,1,2] 


(|Harrison & Kobayashi||2013t 2015 in prep), and ‘r’ and 
‘t' designate reyerse and forward shock, respectively dSari 


fc Piran 

|1999bl 

2003). 

kere we 


•an||1999b[ [Kobayashi fc Zhang||2003[ |/hang et al. 
Here 'we assume that € 3,1 = es^f, where es is 


internal energy density. 


The parameter ^0 describes the shell thickness regime, 
whic h strongly affects the properties of the RS emis - 


sion ( Nakar fc Piran|[2004[|Harrison fc Kobayashi||2013 ). 
Determined as ^0 = where I = 

[(3 —fc)i?/(47rnoi?o^pC^)]^/(^“^^ is the Sedov len gth and 


A n ^ cTqn/G -|- z) is the shell wi dth estimate (Sari fc 


Piran||1995| [Kobayashi et al.||1997|), ^0 ^ 1 and 4o ^ 1 
correspond to the relativistic (thick-shell) and Newtonian 
(thin-shell) regime, respectively. In order to allow for a 
variety of shell thicknesses, we approximate the observed 
^0 distribution using: 


^0 ~ G X (tp.opt/Ygo — 1) 


( 4 ) 


with the value of constant C estimated from hydrody¬ 
namical simulations such that C ~ [5,20,10] for k = 


[0,1,2], respectively (Harrison fc Kobayashi 2013 2015 
in prep). The estimate of ^0 depends on the assumption 
that shell width can be approximated by Tgo; this may 
not be accurate in some cas es (for more discussi on on 
this, see Kobayashi 2000 and Nakar fc Piran|[2004 ). 


optical peaks detected from the literature dRykoff et al. 

2009 

Melandri et al.| 

20101 ILiang et ai. 20131 |Yi et al. 

2013 

Panaitescu et al 

. 2013 Hascoet et al. 

2014p. The 


t To include ^0 in our simulations, we approximated 
e distributions for each k with associated log-normal 
distributions. Although these should not be considered 
as a robust fit due to the sparsity of ^0 values, they re¬ 
semble well the parameter space. Using log-normal dis¬ 
tributions, we generate random values of ^0 to calculate 
the numerical cons tants Cm, Cp, and Ct for k = [0,1, 2] 


(as determined by Harrison fc Kobayashi 20131 2015 in 
prep), which are used in generating the low-trequency 
light curves. 



Figure 1 . Distribution of for ^ = 1 case, as estimated from 25 
GRBs with clearly detected optical peaks (solid black line). Gray 
histogram represents the same distribution, but obtained by taking 
into account also the uncertainty on tp,opt and Tgo- Dashed black 
line is the normal distribution of values, with mean log^o ~ 0.75 
and a ~ 0.42. 


^ Althou gh this assumption m ay not hold for bursts with late 
peaks (e.g. [Guidorzi et al" 1l2014ll. it is still reasonable for many 
observed bursts | |tlascoet et ^ pD14| |. 


Using Equations ([^-(|^, we construct our low- 
frequency model light curves. Light curve scalings in 
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a generalized stra tified med i um us ed b elow are based on 
the derivations of 

time dependencies change 


al. (20141. Both 
according to spectral regime and RS shell type, but we 
initially consider the case where the RS shell crossing has 
already occurred and the scalings are insensitive to RS 
shell thickness. 

2.1. Forward-shock emission 


For a slow-cooling synchrotron spectrum (Sari et 


al. 1998), which is expected at long wavelengths un- 
Her the assumption that r'radio.obs < ^'m.f ~ t'opt at 
ip,opt, the FS peaks at radio frequencies at ip,radio,f 

(^opt/^radio,obs) ^ ip, opt with a peak flux of Fb,niax,f ^ 

(ip,radio,f/ip,opt)“'‘''^^^'‘“''^^i^p,opt, reflecting the e volution 
of Fb.max.f w ith time for cases where fc ^ 0 (e.g. Granot 
& Sari 2002). Here, r'radio,obs and Vopt are the observed 
radio and optical frequencies, respectively. In this spec¬ 
tral regime, the light curve scales as F (x 
and F oc t(-i2p+3fcp-5fc-bi2)/4(4-fc) (Ivan der Horst et al. 


2014), before and after the peak respectively. 


Here we note that a correction factor y = i^m,f/t'opt is 
introduced in cases where Vrn,{ is significantly below the 
optical band, modifying the FS peak tim e by a factor of 
and the FS peak flux by (Melandri et al. 

2010|). Ghanges to RS light curves due to factor y are 


propagated through Equations ([^ and (^. For simplic¬ 
ity, the common assumption to aetect bright afterglow 
emission is y = 1 (i.e. = r'opt), but we discuss 

cha nges to light curves if y = 0.1 or y = 0.01 in Section 

HU 


2.2. Reverse-shock emission 


Using Equations (U and ([^, we can 
peak frequency v^. y and peak flux de 


peaK 


estimate the 
density F^^y, 


of the RS. After the RS shell crossing, Vm,v scales 
as i(i4fc-73)/i2(4-fe)^ causing the peak of the RS light 


mate the RS flux at the observed frequency as: 

(10fc-47)/(12(4-fc)) 




^p,radio,r 


t 


p,opt 


^ . (^radio,obs/^m,r) ^ ^radio,obs ^m,r 

(^radio,obs ^radio,obs ^ ^m,r ; 

(5) 

with the latter half compensating for the possible dif¬ 
ference between the observed radio frequency and the 
peak frequency give n by Equation (§ ( e.g. |Harrison| 
'& Kobayashi 2013). Finally, the light curve in the 


slow-cooling regime scales as F oc t 4/c)/(9(4 fc)) 

F oc t(-73p-bi4fep-b6fc-2i)/24(4-fe) (|van der Horst et al. 
2014), before and after the peak, respectively. 


2.3. Self-absorption approximation 

Synchrotron self-absorption (SSA) becomes impor¬ 
tant at low frequencies. Below the characteristic self¬ 
absorption frequency, Va, which is determined by the 


condition that the optical depth is unity, the system 
becomes optically thi ck and the afterglow flux can be 


reduced significantly (Granot et al. 


1999). The exact 


value of Fa depends on various microphysical parameters, 
which are often vaguely determined. The common alter¬ 
native method used to account for SSA is to approximate 
the SSA flux limit with a simple bl ack body limit with 
the FS or RS effective te mperature (Sari fc Piran||1999b' 


Kobayashi & Sari 2000). It has been sh own that these 
two m ethods do not differ significantly (Shen & Zhang 
2009), and the latter is also more convenient within the 


scope of our paper, as our low-frequency model is based 
on scaling of observed properties rather than using a 
particular model-dependent formulation. At the opti¬ 
cal peak time, the SSA flux limit in the RS region is 
approximated by: 


-^_L 

Abb ^ Q^,geemprQ(l -t- z) [ 


( 6 ) 


where ~ SFoctp^opt is the observed size of the fireball 
and Dl is the luminosity distance. The SSA flux limit in 
the FS region is larger by a factor of Fq due to the higher 
FS temperature. Initially, the limit in the RS (FS) region 
scales as F oc (^p ^2/(4-fc)^ then as 

F oc t(113-22fc)/(24(4- fc)) (p g^\(20-3fc)/(4(4-fc))^ 

Fm. r (Fm f) crosslng (SuH & Piran||1999a| van der Horst 
et al. 2014). We note tiiat the SsA limit will mostly 


(Frr,,.r/FraMn.nhs) k)/(14k 73)^^_^^^^ ( ^ pjj.gg 

• 

1999a 

Kobayashi 

to 

o 

o 

o 

van der Horst et al.||2014). The 

peak 

lux of the J 

iS decays after the RS shell crossing. 

• 


affect the RS emission at early times, but could also affect 
the FS emission at very low frequencies, however not as 
significantly due to a factor of Fq larger limit. 

3. SIMULATIONS 

In order to examine the parameter space of possible 
low-frequency light curves we perform a series of Monte 
Carlo simulations using the model framework of Section 
Each simulation creates 1000 light curves with the 
lollowing input: 

Optical peak time ip,opt: either 200s (‘early’) or 
1000s (‘late’). 

Optical peak magnitude (i? band): either 15mag 
(‘bright’) or 18mag (‘dim’). 

• Observed radio/mm frequency iZradioobs: 1-4 GHz, 
10 GHz or 100 GHz. 

• Circumburst medium structure k: 0, 1, or 2. 

• Circumburst medium density no: lcm“3. We qual¬ 
itatively discuss changes to this parameter in Sec¬ 
tion Hl2l 


• Redshift: drawn from observed normal distribution 
z = 1.84 ± 0.16 ( Salvaterra et al.||2012 |. 


• E, 


74SO- _u_ j _ 

tion logFg^iso = 52.96±0.79 ( Melandri et al.|2014 ). 


drawn from observed log-normal distribu- 


• Shell thickness: Drawn from observed log¬ 

normal distribution (see Figure [B log,Jo ~ 
[0.45,0.75,0.60] ±0.42 (for k = [0,172]). Due to 
sparse distribution, we also tested our simula¬ 
tion using a uniform distribution and the results 
did not change. 
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Micr ophysical parameter s: p ~ 2.36 ±0.59 (normal 


dist., 

dist., 


et al. 


Curran et al.||2010|), Ce € 0.01 — 0.5 (uniform 
Santana et ai.|[20l4 ), ry = 0.2 (fixed, Santana 


2014 r Parameter cb is not explicitly used 


in our formulation because we assume = cs.f 
and use FS properties obtained from optical ob¬ 
servations to normalize FS emission, while the RS 
properties are determined through the shock jump 
conditions (Kobayashi & Zhang||2003 ). 


4. RESULTS 

Figures!^ and [^summarize the light curve predictions 
using the low-frequency model framework detailed in Sec¬ 
tion Figures represent three different observed radio 
frequencies (1.4 GHz, 10 GHz, and 100 GHz). Four pan¬ 
els of each frequency represent simulated light curves at 
three different circumburst medium structures fc, for dif¬ 
ferent combinations of optical peak time and brightness. 
In the following we consider some trends that emerge 
from the collection of light curves. 

In general, it is evident that the need for rapid response 
at low frequencies is extremely important for RS physics, 
similar to optical wavelengths. The brightest and most 
evident RS emission occurs, in most cases, before 0.1 day, 
but for some cases can extend to the order of 1 day after 
the trigger, depending on observed radio frequency and 
circumburst medium structure k. A variety of factors 
contribute to conditions that lead to strong RS emission. 


4.1. Optical peak properties (itp,opt, Tjy.max.f, Vm,t) 


As the deceleration time increases, Fq decreases, caus¬ 
ing the brightness of the RS emission to decrease, but 
also causing the RS and FS peaks to occur later in time. 
Although seemingly counterintuitive, a later deceleration 
time allows for the possibility of catching the decay of 
the RS before it becomes surpassed by the increasing FS 
emission, thus helping to separate the two peaks in time, 
giving each better definition. 

The brightness of the optical peak, and therefore of the 
FS peak, has a rather small overall effect in this model, 
basically just scaling the predicted radio brightness. We 
considered two cases: a ‘bright’ optical peak of 15 magni¬ 
tudes and a ‘dim’ peak of 18 magnitudes. We find that a 
comparatively dim peak, still observable by existing opti¬ 
cal facilities, can make the early-time radio RS emission 
slightly more observable, in the sense that it lowers the 
contribution of the FS emission while still maintaining a 
detectable flux density at low-frequency. 

The optical peak can be caused by the deceleration of 
GRB outflow or by the passage of the typical frequency 
through the optical band. A peak due to the former is 
expected to have a sharp rise, while a peak due to the 
latter would have a slow rise {F cx for fc = 0 and 
shallower for k > 0). If the typical frequency is above 
the optical band at the deceleration time, the optical 


flux is expected to rise as t^^^ or shallower until the typ¬ 
ical frequency passes through the optical band. Since 
such a s low rise is hardly observed (e.g. Melandri et 
al. 2010|, we assume that the typical frequency is 


already below the optical band z/opt at the deceleration 
time: x < 1. Ghang^ when x is signihcantly below 1 are 
presented in Figure H (left), where we show the evolution 
of light curves for a case of fc = 0, tp,opt = 1000 s and 


R = 18 mag. The peak times of both FS and RS emis¬ 
sion are shifted to earlier times, while peak fluxes are 
shifted to higher values, but especially in the RS regime 
flux density evolution is masked by the SSA limit, for 
which the break also shifts due to shift. Lower y 
thus makes RS emission less pronounced. 


4.2. Circumburst medium 

Different values of circumburst medium stratihcation 
k affect the predicted light curves similarly to decreasiM 
the optical peak flux. As shown in Figures and 
the FS peak becomes less pronounced as we move from 
fc = 0 to A: = 2, due to the —fc/(2(4 — k)) dependence of 
the peak flux. The slope before the FS peak, however, 
becomes more shallow as k increases, contaminating the 
early-time RS emission. 

Gircumburst medium density ug determines the 
Lorentz factor estimate (Equation and consequently 
affects RS emission properties (EquatioiisJ^ and and 
the SSA limit approximation (Equation pT Due to the 


eling of afterglows 

(Yost et al. 

20031 Cenko et al.|2011 

Japelj et al. 2014) 

which cou 

d span over 6 orders ol 


our simulations as otherwise the spread of the simulated 
light curves would be too wide to identify key trends 
and effects. However, we do assess the evolution of light 


curves when changing ng from 0.001cm ^ to 1000 cm 

as obtained from 

Genko et ah 

(2011). The resulting light 

curves are plotted in i:^‘igure|4 

(right), tor a case ot fc = 0, 


Ap.opt = 1000s and R = Isrnag. We infer that in this 
case, the RS peak is much more pronounced in envi¬ 
ronments with lower circumburst medium density, while 
shifted towards earlier times. Same results also hold for 
different values of circumburst medium structure k. 


4.3. Self-absorption and zZradio,obs 

The observed radio frequency has the two-fold effect of 
normalizing the level of the SSA flux and influencing the 
peak time of both RS and FS. The effect is more apparent 
in the latter as it is dependent on the time evolution 
of while the RS will peak comparatively early due 
to its Fg dependence. Around ^ 1 GHz, SSA tends to 
dominate and severely reduces the RS flux level, in some 
cases even removing any temporal structure that would 
help identify the RS emission. Since the SSA flux limit 
is oc zz^, the effect is reduced with increasing observing 
frequency, aiding the observation of the RS at > 10 GHz. 
However, when deciding on an observing strategy, it is 
necessary to balance these opposing effects by choosing a 
frequency that will not be too self-absorbed nor too high; 
in the latter case the RS and FS will peak very early and 
the RS signature would be very difficult to detect. 

4.4. Numerieal corrections 

Further folded into our simulations are the changes in 
the light curves due to the addition of the numerical cor¬ 
rection factors in Equations 0 - Particularly at 

early time, many of the figures showoands of predicted 
light curves covering several orders of magnitude in flux 
density. This effect originates partly from the wide dis¬ 
tribution of Ce and its effect on the SSA limit but also 
from the numerical constants. We introduced these fac¬ 
tors by assigning each randomly generated light curve a 
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t [days] 


t [days] 


Figure 2. Predicted light curves for 1.4 GHz. Thick lines represent the median light curves, while thin lines represent simulated light 
curves (Section!^ we randomly chose 100 out of 1000 simulated light curves for clarity). Colors indicate different circumburst medium 
stratification /c. i^our panels represent various peak times tp,opt and peak magnitudes R. Also shown for /c = 0 cases are the relative 
contribution of the RS corrected for corresponding SSA limit (dashed black line), and of the FS corrected for corresponding SSA limit 
(dotted black line). 


value for the dimensionless RS shell thicknes^arameter, 
^ 0 , from the observed distributions (Figure [l for /c = 1 
case). The effects in intermediate (^o 1) and thin-shell 

(^0 > 1) regions are particularly pronounced a nd on the 


scale of 2 orders of magnitude (see Figure 3 in Harrison 


& Kobayashi 2013 for fc = 0 case). In practice, this ef¬ 


fect can be mitigated since can be approximated from 
the deceleration time and burst duration (Equation . 
This approximation, however, assumes that the emitted 
shell has a homogeneous structure and that the deceler¬ 
ation time occurs after the end of the prompt emission. 
These are reasonable assumption s but may not alway s 
be valid, such as in GRB 061126 (Gomboc et al. 2008). 
This adds further uncertainty to tne low-^o value tail of 
the distribution. 

The broad distribution causes a large spread of RS 
light curves but also provides parameter spaces in which 
the RS can more easily be observed. We can better un¬ 
derstand the effect that these numerical factors have by 
looking at the values of burst parameters that produce 
bright (dim), or above (below) median flux light curves 
in our simulations. Figure!^ shows the and E^^iso dis¬ 
tributions for an example trial of 10 GHz light curves 
from an assumed 18**' magnitude optical peak occurring 
at 1000 s after GRB. The brighter light curves are created 
by a combination of higher E^^iso and smaller ^o^ since 
the former contributes to a larger Fq and the latter’s 


effect is less pronounced closer to the thick-shell regime. 

5. DISCUSSION 

Our simulations show that later optical peaks, assum¬ 
ing they indicate the hreball deceleration time, frequen¬ 
cies between > 1 and < 100 GHz, and bursts that have 
a circumburst stratification closer to fc = 0, are preferred 
candidates for observing RS emission. Brightness of the 
optical peak has less of an effect in this model, although 
this is also dependent on the burst parameters and ob¬ 
serving frequency, as the RS component may be partially 
self-absorbed or mixed with emission from the FS. RS 
emission with median flux levels of the order of /rJy- 
mJys should be visible at early times (< 0.1 — Iday) 
with respect to the GRB trigger and are observable with 
currently operating facilities. Future facilities will en¬ 
able shorter exposures times due to improved sensitivity, 
allowing for better light curve sampling. 

As an example, with a standard array configuration 
(34 1-m, 9 7-m, and 2 TP antennas) observing at 100 
GHz (ALMA Band 3), ALMA can detect afterglow emis¬ 
sion (3cr) to limits of 10 (1, 0.1, 0.0^ mJy with 0.06 s 
(5.84 s, 9.7 min, 37.4 min) integration^ Modest integra¬ 
tion times can lead to the detection of many of the bright- 

^ Calculated using the ALMA sensitivity calculator found at 
http: / / almascience.eso.org/proposing/sensitivity-calculator 



























F [mJy] F [mJy] F [mJy] F [mJy] 


6 


Kopac et al. 



10 ^ 

t [days] 


10 ^ 

t [days] 
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t [days] 


Figure 3. Predicated light curves for 10 GHz (top 4 panels) and 100 GHz (bottom 4 panels). See description in Figure]^ 
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t [days] 


tp.opt = 1000 s; R = 18 mag 


t [days] 


Figure 4. Evolution of light curves (fc = 0, tp,opt = 1000 s, R = 18 mag case) at 10 GHz (gray) and 1.4 GHz (magenta) when changing: 
(left) X = f/^'opt parameter (x = 1, 0.1, or 0.01, from bottom to top). Smaller x shifts FS and RS peaks towards earlier times, while 
peak fluxes are shifted to higher values but are masked by the SSA limit especially for the RS emission (dashed lines). Simulated light 
curves (thin lines) are only plotted for x = 0.1 case, (right) Circumburst medium density (no = 0.001, 0.1, 10, or 1000 cm“^, from left to 
right corresponding to RS peaks; for n = lcm“® case see Figures and |^. Larger no lowers Fq estimate, shifting the RS peak towards 
later times and lower flux densities (dashed lines). Simulated light curves (thin lines) are only plotted for no = 0.1 cm“^ case. 



Figure 5. Distribution of ^o (left) and (right) for a trial simulation (tp,opt = 1000s, peak mag = 18mag, observed frequency = 

10 GHz, fc = 0) for fluxes above and below the median. Bursts with above-median fluxes in the radio have larger energies and smaller ^o 
(closer to thick-shell regime), helping to mitigate the effects of the numerical corrections. 


est afterglows, and in the k = 0 and k = 1 cases, a larger 
fraction of the dimmer afterglows. Longer wavelengths 
that are more affected by self-absorption will pose big¬ 
ger problems, but should be observable in the extreme 
cases. In addition to sensitivity, facility response time 
is extremely important and capturing clear RS emission 
will require aggressive ToO campaigns. 

Polarization can provide further evidence of RS emis¬ 
sion. This technique has been pioneered at optical wave¬ 
lengths and is rapidly becoming feasible at radio and 
sub-mm wavelengths. Polarization of GRB emission is 
a prediction of various theoretical models and has been 


observed in prompt y-rav emission dGotz et al. 12009 

Yo- 

netoku et al.||20I2 

Gotz et al. 112014) and in optical 

to a 


level of 10-30% in the first minutes after the burst (Steele 
et al. 12009 Uehara et al.|2012 Mundell et al.|2013 1, prob¬ 
ing the tundamental fireball magnetic held. 


magnetic 

At later times (hours to days post burst), low levels 
of optical polar ization of a few % have been reported by 
various works (IGovino et al.||1999| IGreiner et al. 1 12003 


Wiersema et al^L)12[ 2U14[ ) and are thought to be due 


to h'S-dominated emission, probing the shocked ISM. Ra¬ 
dio polarization for a GRB has yet to be detected and 
upper limits of a f ew percent have been reported for a 


handful of b ursts (Taylor et al. 


2004 2005 Granot & 


Taylor 20051 at days alter the GrB trigger. Although 


emission at such late times most often comes from the 
FS, the underlying component of RS emission would help 
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us understand its effect on the overall ligh t curve behav¬ 
ior and pola rization properties. Recently, |van der Horst, 
et al. (|2014 1 undertook radio polarization observations of 
GRB 130427A at 1.5 and 2.5 days, finding similar upper 
limits despite the fact that these observations were taken 
at the radio peak time, which was associated with the RS 
emission. This event further stresses the importance of 
rapid multi-wavelength observations if we are to under¬ 
stand the size and scale of the magnetic fields in GRB 
ejecta. 

Looking beyond current capabilities, the SKA will pro¬ 
vide unpreced ented sky coverage at frequencies around 
1.4 GHz (e.g., Ghirlanda et al.|[2013 Burlon et al.||2015 l 
and has the potential to observe the very early-time RS 
signature of GRB emission. However, at such low fre¬ 
quencies, the early time flux in many cases will mostly 
be suppressed by self-absorption (see Figure!^ and ex¬ 
pected RS flux densities at 0.1 days after the^RB trig¬ 
ger are on average lower than < 10 /iJy. But depending 
on burst’s parameters, it is still possible in some cases 
to detect distinguishable RS peaks with peak flux up to 
^0.1 mjy, as indicated especially by thin lines from sim¬ 
ulated light curves (especially for ‘dim’ and ‘late’ optical 
peaks occurring in k = 0 type medium, see Figure]^. 
Such fluxes could be achieved by the SKA with a reason¬ 
able (^ 10 — 30 minutes) integratio n times with a typical 
array setup (Ghirlanda et al. 2013). This exposure time 
would be short enough to achieve a reasonably deep limit 
without compromising valuable temporal resolution. 


6. CONCLUSIONS 

We have presented a collection of predicted radio after¬ 
glow light curves in the framework of the GRB standard 
model, considering contributions from both FS and RS 
emissi on within an updated low-frequency model frame¬ 
work (Mundell et al. 2007 Melandri et al. 2010). Our 
results show that: 


• Inclusion of correction factors from numerical sim- 
ulations that parametrize sh ell-thickness regime 
(Harrison & Kobayashi 2013 2015 in prep) pro- 
vides more accurate estimates of RS emission. 


• RS emission can significantly alter the temporal be¬ 
havior at e arly times, producing a distinct “radio 
flare” (e.g. |Kulkarni et al. 1999). The brightness 
and shape depend on various parameters, but are 
mostly affected by the SSA, which is stronger at 
lower frequencies. 


• RS emission could be best identified at frequen¬ 
cies around ~ 10 GHz, peaking at < 0.1 day af¬ 
ter the burst, depending on burst’s parameters. 
At frequencies around ~ 1 GHz, RS emission is 
more likely to be strongly suppressed by the SSA, 
while at higher frequencies (around ^ 100 GHz), 
RS emission tends to peak at earlier times (< 0.01 
day after the burst). 

• Gurrent and future radio facilities with high sen¬ 
sitivity and short response times could detect ra¬ 
dio flares with peak flux densities around F ~ 
0.01 — 0.1 mJy {F ^ 0.1 — lOmJy) at t < 1 day 
(t < 0.1 day) after the burst, at 1.4 GHz (10 GHz), 
respectively. 
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